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Abstract 

This paper gives sufficient conditions for having complete synchroniza- 
tion of oscillators in connected undirected networks. The considered oscil- 
lators are not necessarily identical and the synchronization terms can be 
nonlinear. An important problem about oscillators networks is to deter- 
mine conditions for having complete synchronization that is the stability 
of the synchronous state. The synchronization study requires to take into 
account the graph topology. In this paper, we extend some results to 
non linear cases and we give an existence condition of trajectories. Suffi- 
cient conditions given in this paper are based on the study of a Lyapunov 
function and the use of a pseudometric which enables us to link network 
dynamics and graph theory. Applications of these results are presented. 
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1 Introduction 

The study of the dynamies of coupled nonlinear dynamical systems are the 
subject of a growing interest in various communities like in theoretical physic, 
in information technology or in neuronal biology. The literature on this topic 
shows different kinds of synchronization (see [10]). Classically, two coupled 
limit-cycle are said synchronized when their time evolution is periodic with the 
same period and perhaps the same phase. From the discover of synchronization 
of chaotic systems (see [1, 5, 8]), the word synchronization recovered different 
meanings such as having identical or functional related solutions, eventually with 
a delay. The definition has also been modulated by considering strong forms 
like complete, cluster form or weaker forms like phase and lag synchronization 
(see [11]). 

An important question about synchronization of a network of oscillators is 
to determine the stability of the synchronisation state. This question leads to 
consider some properties of networks and state vectors of oscillators (see, for 
example, [4, 13, 14, 15, 17]). For this purpose, two methods are proposed in 
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the literature. The first one called master stability function is based on the 
computation of a Lyapunov exponent and the eigenvalues of the connectivity 
matrix [:;]. However, this method is adapted when the coupling terms are lin- 
ear and the computation of eigenvalues can become a difficult task. A second 
proposed method is the connection graph stability method (see [4]). It links the 
study of a Lyapunov function and the graph topology. This productive method 
has been extended to unbalance and undirected graph (see [2, 3]). 

The results presented in this paper generalize some results of [4] to the non 
linear synchronization case. For this, we introduce a notion of pseudometric in 
the graph. The determination of the sign of the Lyapunov function derivative 
requires two steps. The first one is to use assumptions allowing comparisons 
between oscillators and synchronization terms. The second step consists in 
using pseudometrics which enable us to use some graph properties. For the 
complete synchronization, we present two results. The first one gives a condition 
on synchronization strength for having a global synchronization of oscillators. 
The second result is a local versus of the first one, that is when the oscillators 
are closed to the synchronization variety. In these two cases, we give sufficient 
conditions that insure existence of trajectories. 

This paper is organized as follows. The problem statements are presented in 
Section 2. First, we precise the kind of systems and the kind of synchronizations 
considered. Then, we recall the definition and some properties of pseudomet- 
rics defined on a graph. In Section 3, after precising the assumptions on the 
synchronization term, main results, that is conditions for having complete syn- 
chronization of the system of oscillators, are presented. These results are applied 
in Section 4. 



2 Problem statements 

Thereafter, is the transpose of the vector Y = {Y^,..., F™) € R". 

2.1 Systems and synchronizations considered 

Let G be a connected undirected graph and n its number of vertex. The graph 
G describes the set of interactions between the oscillators. We denote by £ the 
set of its edges. If G contains an undirected edge from a vertex z to a vertex j, 
we denote it by (i, j). 

The considered dynamical systems are defined by the following system of 
equations: 

' Xi = Fi(Xi,t) - e ^ h{Xi,Xj), 

(lj)6£ 

; (1) 

Xn = Fn(Xn,t) — e h{Xn,Xj), 
{n,3)e£ 

where 

• Xi = {Xl, . . . ,Xf)'^ is the vector composed of the d coordinates of the 
i-th oscillator, 

• Fi = (Fl, . . . , Ff)'^ is the vectorial function defining one oscillator, 
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• h = {h^, . . . , h'^Y' is the synchronization function which defines the vector 
couphng between osciUators, 

• the real parameter e corresponds to the synchronization strength 

RecaU that, for a given initial state of the set of oscillators (Xi(0), ^2(0), • • • Xn{0))'^ , 
system (1) synchronizes completely if, for all G 

|jX,(t)-X,(t)|| — — ^0. 

This means that the vector (Xi, . . . , X„) approaches the synchronization man- 
ifold defined by Xi{t) = X2{t) = • • • = X„(t). In particular, this implies that 
the oscillators have the same asymptotic behavior (such as chaotic trajectories, 
stable and periodic solutions). The complete synchronization of all oscillators 
can occur whatever their initial states are, in this case, the synchronization is 
said global; otherwise it is said local. 

In this paper, we focus naturally on the differences A.^^ = Xj — Xj and 
therefore on the vector 

A = (Ai,2, • • • , Ai,„, A2,3, • • • , A2,„, • • • , A„_i,„)^ . 

Thus, proving the complete synchronization of system (1) is equivalent to prove 
that ||A(0|| ^0. 

2.2 Quasimetrics defined on a graph 

In the following, we consider pseudometric verifying the p-relaxed triangle in- 
equality for a positive real p, that is an application (p : D x D ^ R+, where D 
is an non empty set, satisfying the following three axioms: 

• ip{zi,zi) = 0; 

• ip{zi,Z2) = ip{z2,zi) (symmetry property); 

• Lp{zi, Z3) < p {(p{zi, Z2) + f{z2, z^)) (p-relaxed triangle inequality). 
Remark that any classical metric is such a pseudometric with p = 1 . 

Let iy9 be a pseudometric on a set D. Let's set, for all m G N*, p{m) the 
smallest real such that 

ip{zi,Zm+l) < p{m) [ip{zi,Z2)A Lp{Zyn, Zra+l)] ■ (2) 

Note that p{l) = 1. 

In the following examples, expressions of p{m) appearing in inequalities (2) 
are direct consequences of the convexity of functions x — > (a;^)" and x — )■ 

Example 2.1. 1. The application (yS^ : x K+ defined by 




with a > 1/2 is a pseudometric for which p{m) = rn^' 
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2. Let D he the closed hall of center and radius 2 — a/2- The application 
ip : D X D defined hy 



Xi \ 




( X2 




M 


2/2 









is a pseudometric for which p{m) = m. 
We have the foUowing properties. 
Proposition 2.1. 1. The sequence of reals {p{m))m>i is increasing. 

2. For all m e W , we have p{m) < p™^^ (see [16]). 

3. Let ifi and ip2 he two pseudometrics on D and pi{m) and P2{'>n) he 
the smallest respective reals verifying (2). For all a > and /3 > 0, 
the application a ifi + (3 ip2 is a pseudometric on D satisfying p{m) — 
Max{pi (to) ,p2{m)}. 

We now apply pseudometrics to networks of oscillators. Recall that a state 
vector Zi of an oscillator is associated to j-th vertex of G. Let's consider a 
pseudometric Lp defined on the set of state vectors of oscillators. This pseu- 
dometric enables one to define the pseudolength ip(zi,Zj) between vertices i 
and j and also the pseudolength ip{zi-^ , Zi.-^) + • • • + ip{zi^_-^, Zi^) of any path 
Pi,j — (i = ill *2, ■ ■ • ,im = j) from vertex i to vertex j . 

In the following proposition, we bound, up to a multiplicative constant C(G), 
the sum of pseudolengths between any two oscillators by the sum of pseu- 
dolengths of paths joining any two oscillators. This constant plays an important 
role in Theorems 3.1 and 3.2 since the synchronization strenght e appearing in 
these theorems is proportionnal to this constant. 

Proposition 2.2. Let G he a connected graph, £ he the set of its edges and ip he 
a pseudometric on a set D . For any vertex i, let Zi € D he a vector associated 
to vertex i. There exists a constant C depending only on G so that we have 

^ip{z,,Zj) <C ^ 'p{z^,ZJ). (3) 

Moreover, the smallest real C satisfying (3), C{G), is hounded hy 

"^^^mpim), (4) 

where 5{G) is the diameter of G. 

Proof. Let i and j be two vertices of G and let's denote 
Pi.j = {i = ii,i2, - ■ ■ ,is+i = .j) 

a path of G from the vertex i to vertex j (recall that G is connected). Since ip 
is a pseudometric on D, we have <p{zi, zj) < p{s) J2e=i 'fii^in ^ii+i) ■ 

The path Pi^j can be chosen so that s < d{G). Suppose that this choice is 
done for any vertices i and j; since the sequence {p{n))neN* is increasing, we 
have p(s) < p{S{G)). Consequently, for any vertices i and j, we have (p{zi, zj) < 
p{S{G)) 5{G) Max {{(p{zi, Zj) \ {i,j) G £}) which implies the result. □ 
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In Theorem 3.1, we need to determine the lowest bound C(G) of the set of 
reals C satisfying inequality (3). The bound (4) of C{G) may not lead to a good 
estimation of C{G) for a particular graph; nevertheless, this bound is valid for 
any graph with n vertices. 

In the case of a pscudomctric satisfying the classical triangle inequality, i.e. 
when p{n) — n for all n G N*, a method taking G as input and returning a 
bound of C{G) is proposed in [ !]. Its two main steps are: 

1. for all (i, j) with i > j, choose a path Pi^j', this path is usually chosen with 
minimal length (number of edges in the path) ; 

2. for each edge e of the connection graph, determine the sum B{e) of the 
lengths of all chosen paths Pij containing e. A bound for C(G) is then 
Max{B{e) : e e £}. 

For each choice of paths, these two steps return a bound for C{G). Clearly, 
the number of possible paths is huge but computations of bounds for C(G) are 
possible since most of these choices are suboptimal. Up to a slight modification 
of the first step, this method can be applied here: its consists in considering, 
for all path P^j-, the pseudolength p{\Pi,j\) instead of its length \Pi.j\- 

Remark 2.1. In the case of pseudometrics ip satisfying p{m) = m, explicit 
hounds of C{G) for specific graphs and the method proposed in [,{, :i] for com- 
puting C{G) from G can be directly used. This is the case of the second function 
in Example 2.1. 

3 Complete synchronizations 
3.1 Hypothesis 

Afterwards, two cases are considered. The first one is the global complete 
synchronization for which oscillators Xi,. . . ,X„ lies in D = W^. The second 
one is the complete synchronization for which oscillators are in a neighborhood 
D of the variety Xi = X2 = • • • = Xn ■ 

Thereafter, we will suppose the following assumptions on system (1). 

• For all G £, there exist some non negative reals ai, . . . , Od such that 

d 

V(X„ X,) e D, X,) - """{X^ - Xf)h\X,, Xj) (5) 

k=l 

are pseudometrics where h — {h^ , . . . , h'^)'^ is the synchronization func- 
tion. 

• For all £ |1, n]^ and, for all t > to where to e R, 

d 

V(X„X,) e D, ^a,(Xf - X^) {Ft{X,,t) - F^{X,,t)) < v{X.„X,) . 
fc=i 

(6) 
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. For all e II, nf, V(X„X,) € D, 

^{X,,X,) = and/or ELi ^ (^"'(^Mt) - = 



Remark 3.1. 1. Notice that hypothesis (5) implies that, 

y(i,j) G £, \f{Xi,Xj) e D, h{Xi,Xj) — -~'h{Xj, Xi) (antisymmetry). 

(8) 

2. The assumption (7) is necessary for proving the complete synchronisation 
of system (1) in Theorems 3.1 and 3.2. The condition ip{Xi,Xj) = m 
this assumption is not always sufficient when it does not imply equalities 
of all the components of oscillators. In this case, the second condition is 
necessary for proving the complete synchronization. 

For practical cases, a first problem is to prove the existence of trajectories of 
system (1) for a sufficient large t. For this goal, the following proposition enables 
us to link existence of trajectories between synchronized and non synchronized 
systems. 

Proposition 3.1. For all € [l,nf, suppose that assumptions (5), (6) 

and (7) are satisfied and that, for all t > to, 

Xj'F,{X,,t)<^{\\X,\\) 

where ^' satifies the conditions 
/■+°° ds 

/ —— = +00 and ^(s) > for all s > sq > 0. 

J s=so ^V') 

Then, the Cauchy's problem defined by system (1) and an initial condition 
Xiito) \ 

• g j^nd f^^g ^ solution on the complete semi-axis [to;+oo) . 

Xnito) J 

Proof Let's set X= : e R'"* and = : eK'"^. In 

V ^„ / V F,,{Xn,t) j 

a first step, we prove that there exists a real /3 such that the following inequality 
between the scalar products holds: 

X'^X < px'^F{X,t). (9) 

For this, we consider the dnxdn diagonal matrix M = Diag{ai, . . . a^, ■ • ■ , oi, . . . ad). 
We have: 

n d n d 

X^MX = ^^afeX,'=i^'=(X„t)-e^^afe ^ X^h'^{X,,X,) 

1=1 fc=l 1=1 fc=i {j\(,i-j)e£} 

d 

= X^MF{X,t)^eY, akX'^h''{X,,X,) 
fe=i (ij)e£ 
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and, since to any edge € £ corresponds the edge (j, i) S we obtain 

d 

XT MX = X^AfF(X,t)-|^afc ^ X^h'^iX^^Xj) + X'^h^iX.^Xi) 

fe=i (ij)e£ 

d 

= X'^MFiX, - I ^ afc ^ - X,) (see equality (8)) 

k=l (ij)6£ 

< X'TMF{X,t). (see assumption (5)) 

Inequality (9) is then a direct consequence of the fact that the reals are non 
negative. 

If the conditions of the proposition are verified, inequality (9) shows that we 
have, for alH > to, _ 

x^i:<vi/(||x||) 

where 4* is a application satifying the conditions 
/■+°° ds 

/ - — = +00 and *(s) > for all s > sn > 0. 

Thus, system (1) satisfies the conditions of Wintner's theorem ([12]) and, 
consequently, solutions of system (1) are defined for any t > to- 

□ 

3.2 Global synchronization 

Theorem 3.1. Suppose that the assumptions done in Section 3.1 are satisfied 
for D = (R'')^. If e > — — , where Cq is the optimal bound such that inequa- 
lity (3) holds, then system (1) synchronizes completely. 

Proof. In order to show this result, we will apply the second method of Lya- 
punov. Let's consider the Lyapunov candidate function: 

k=l i<j 

Clearly, this function is non negative if A 7^ and equal to iff A = it that 
is when the system (1) is synchronized. 
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The derivative of V gives: 
V 





d n 71 
fc=l i=l j = l 


d 

fc=i 


' n n 


d 

fc=i 


n[j2x':Fi^{X,,t)~ 



»=1 {il(«j)G£} 



E«fc 

fc=l 



E«fc 

fc=l 



/ (»j)e£- 



Since each edge e f corresponds to an edge (j, i) and using equahty (8), 
we have, for ah fee [1, Jt-], 



2 Y h\X,,X^) 



^ h\X,,X,)+ Y -h\X,,Xj) 





and 



2Yak Y X^h^{X,,X,) = Y 



ak 



k=l 
d 

E^fc 



fc=l 



^ x1h'^{x,,x,)+ Y xfh^ix,,x,) 



Y X^h\X,,X,)+ Y 'Xfh\X,,X,) 

^ (^(X„X,)(see5). 
(ij)e£ 
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Moreover, we have 

These three equahties gives 

d 

V = J2Y.Y'^^'-^'HP^HX^,t)^F^{X,,t))-ne ^ ^{X,,X,) (10) 

i,3 k=l {i,j)eS 

With assumption (6) and incquahty (3), we obtain 



Since ip is a pseudometric the right factor of this last expression is non negative. 
Cg 

Therefore, if e > then V < 0. To prove that V is negative definite, it 

2n 

remains to show that if y = then Xi = X2 = • • • = X„. Suppose that V = 0. 
Since — ne) < 0, the last inequality implies that we have ip{Xi,Xj) = for 
all G £■ From equality (10), we obtain 

d 

i,j k=l 

Consequently, assumption (7) is satisfied and system (1) synchronizes. □ 
3.3 Local synchronization 



Let H be the diagonal matrix Diag{ai, . . . , ad) and H 



the matrix composed with matrices H. The application 



/ H ••• \ 
F ••• 



V ■■■ H J 
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X ^IXT-RX 
is a norm since ai, . . . , are non negative. Let's set 

v{t)^\\mrv = \jz E «^(xf(i)-^i=wr 

A;— 1 i<j<7i 



(11) 



9 



Theorem 3.2. Let B the closed hall {X G W^^^'^ \ \\X\\v < r} where r is a 
non negative real. Suppose that assumptions of Section 3. 1 are satisfied when 

o o 

A belongs to the inner B of B and suppose that, for an instant tg, A(to) G B- 
Cg 

If € > — — , where Cg is the optimal bound such that inequality (3) holds, then 
2n 

system (1) synchronizes. 

Proof. Let's show that if A(io) &B then Vi > t^, A{t) e B. If A(to) ^B, by 
definition of B, we have V{to) < r^. Suppose that there exists ti > to such that 
A{ti) ^ B; by definition of B, we have V{ti) > r^. Since t V{t) is continuous, 
there exists a real t2 = Inf{t G [to,ti]\V{t) = r^}. The mean value theorem 
shows that there exists ^3 G (^0)^2) such that V'its) = > 0. 

On the other side, since <3 < t2 = Inf{t e [to,ti]|y(t) = r^}, we have V{t3) < 

o 

r and A(t3) G B- Consequently, the hypothesis of Section 3.1 are satisfied 
by A (is) and wc can proceed like in the proof of Theorem 3.1 to show that 
^'(^3) ^ 0. This brings to a contradiction. 

Finally, we have Vi > to, A(t) G S and the assumptions of Section 3.1 are 
satisfied for any t > to. Now, we can proceed like in the proof of Theorem 3.1 
to conclude. □ 



4 Applications 

In this section, we focus on applications of Theorems 3.1 and 3.2 in order to 
have a sufficient condition for global synchronization of two systems. The fact 
that solutions of these two systems are defined on M is a direct consequence of 
Proposition 3.1. 



4.1 Global synchronization of a network of neurons 

In this section, we apply Theorem 3.1 to a network of neurons satisfying the 
FitzHugh-Nagumo model (See [G]). Recall that the dynamic of a single neuron 
is modelised by the equation X = F{X) where 



• X = 



some real parameters a, b, c and d. 

In the following, we suppose that b is positive. Let's set G the connected 
graph describing the interaction between the oscillators, n its number of vertices 
and £ the set of its edges. For the synchronization terms, we consider the 
function h defined by 



V(i,j) e 11, nf, h{X„Xj) = '''^{y-V)' 
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with a > 1, /3 > and 7 > A4ax{0,—c}. The system of equations for the 
network of oscillators is then 

' li=Fi(Xi)-e ^ h{Xi,Xj), 

: (12) 

Xn — Fn{Xn) i h{Xn,Xj). 

(n,j)e£ 

The three hypothesis of Section 3.1 are satisfied with ai = 1 and 02 = 1/6. 
Indeed, 

1. assumption (7) is obvious; 

2. the fact that the application ip corresponding to h, explicitly defined by 

ip{X,,Xj) = a{xi - Xj)^ + P^{x,- Xjf +j/b{y, - y^)^ , 

is a pseudometric satisfying p(m) = m^^^ is a consequence of Example 2.1 
and Proposition 2.1. Therefore, assumption (5) is satisfied; 

3. the following inequalities shows assumption (6), for all {Xi,Xj) e D, 

ELi - X^) (FtiXi) - F^{X,)) 

Xi - Xj \ / -(x^i - x^j) + {xi - Xj) - {yi - yj) 
Vi -Vj J \ {xi -Xj)- c/b[yi - yj) 

= - Xj){x^ - x]) + {xr - Xjf - c/b{y, - yjf 
<V{X,,X,). 

For any connected graph G with n vertex, inequahty (3) is verified for the 
bound of C(G) given by C = — -S{G) p{d{G)). Theorem 3.1 shows then 

(n — 1) (5(G)®/'^ 

that, for any connected graph G with n vertex, if e > ^ then 

system (12) synchronizes. 

4.2 Local synchronization of a network of oscillators 

In this section, we apply Theorem 3.2 to a network of Chua oscillators. We 
consider the simplified version suggested by Chua for these oscillators (see [7]): 
if we set X ~ {x,y,z)'^, the state equation for a single oscillator is given by 
X = F{X) where 

a[y-x- f{x) 
F{x,y,z) = \ x-y + z 
—by — cz 

a>Q, b>0, c>Q and / is a piece- wise function f{x) = dx + l/2{d — e)(|x + 
1| - |a; - 1|) with 2d < e. 

Since / is a piece-wise function, a real S > bounds the set of slopes 
|/M^/(y) I < |a; — y| < l|. In the following, we suppose that: 
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1. the set of vertex of G is f = {(1; 2), (1; 3), . . . , (1; n)}. In other words, we 
consider a star configuration of oscillators; 

2. the synchronization function h is given by 

h{{xi,yi,Zi),{xj,yj,Zj)) = \ 



The equation for the i-th oscillator of the network is then 

\ ( a[yi - -^i - fix-i)] \ / ad{xi - Xj)e^^^'''^''^ 

2/i = x,-y^ + Zi + *^ X! ^ 

Zi J \ -byi - czi J j I {i,j)e£ \ 

Assumptions of Section 3.1 have to be verified in order to apply Theorem 3.2. 
The first one is obvious. For the second and the third one, let's set ai = 1/a, 
a2 — 1 and 03 = 1/6. 

Let's consider a closed ball 6 = |x e R^^^'^ \ \\X\\v < (V^ - l)Vaj 
where ||.||y is defined by (11) and the norm given by 



where H is the diagonal matrix Diag{ai, . . . ,ad). If we have A € i3 then 
II Aijily < (v^ — l)-y/a. This implies that \ Xi — Xj |< 2 — and, according to 
Example 2.1, the application (f corresponding to h satisfies assumption (5). 

Let's verify assumption (6). We have 

/ Xj-Xj \ 



a 

Vi - Vj 



4iy^ - yj) - {^^ - - (fix,) - f{xj))] 

(x, - Xj) - (y, - yj) + {z, - Zj) 
-KVi -Vj)- c{z, - Zj) 



= {x^ - Xj){f{xi) - f{xj)) - {xi - Xj)'^ - {yi - yj)^ - c/b{zi - Zj)' . 
By definition of S, we have {xi — Xj){f{xi) — f{xj)) < S{xi — Xj)'^e^^^^^~^^^ . 
This shows inequality (6). 

Moreover, if ^{x,,x,) = and ELi ^kiX^ - X]) {F^{X,) ~ F^{X,)) = 
then we have Xi = Xj. Consequently, assumption (7) holds. 

Since the induced pseudometric satisfies Vm G N*, p[m) = m (see Exam- 
ple 2.1), the bound Cq is given explicitly by 2n — 3 (See Remark 2.1 and [4]). 

o 

Theorem 3.2 can now be applied : if A(to) € B for an instant to and if 

2n-3 , , , , . 

e > then system (1) synchronizes. 

2n 
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5 Conclusion 



In this paper, sufficient conditions for proving complete synclironization of os- 
cillators in a connected undirected network are presented. The contribution of 
this paper lies in the extension of results established in the case of linear syn- 
chronization to the non linear case. For this, we have introduced pseudometrics 
which enable us to link graph topology and minimal synchronization strength 
between oscillators. Under our assumptions, a criterion proving the existence of 
trajectories is given. Two results for proving the complete synchronization are 
then proposed: the first one gives a global criterion and the second one deals 
with local synchronization, that is when the trajectories lie in a neighborhood 
of the synchronization variety. To illustrate these results, two applications are 
treated. 
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